Unconventional phase transition of phase-change-memory materials for optical data storage
Chen Nian-Ke, Li Xian-Bin
State Key Laboratory of Integrated Optoelectronics, College of Electronic Science and Engineering, Jilin University, Changchun 130012, China

 

† Corresponding author. E-mail: lixianbin@jlu.edu.cn

Abstract

Recent years, optically controlled phase-change memory draws intensive attention owing to some advanced applications including integrated all-optical nonvolatile memory, in-memory computing, and neuromorphic computing. The light-induced phase transition is the key for this technology. Traditional understanding on the role of light is the heating effect. Generally, the RESET operation of phase-change memory is believed to be a melt-quenching-amorphization process. However, some recent experimental and theoretical investigations have revealed that ultrafast laser can manipulate the structures of phase-change materials by non-thermal effects and induces unconventional phase transitions including solid-to-solid amorphization and order-to-order phase transitions. Compared with the conventional thermal amorphization, these transitions have potential superiors such as faster speed, better endurance, and low power consumption. This article summarizes some recent progress of experimental observations and theoretical analyses on these unconventional phase transitions. The discussions mainly focus on the physical mechanism at atomic scale to provide guidance to control the phase transitions for optical storage. Outlook on some possible applications of the non-thermal phase transition is also presented to develop new types of devices.

1. Introduction

Phase-change memory (PCM) was proposed by Stanford Ovshinsky in 1968.[1] The information storage of PCM relies on the big contrast of the electronic or optical properties between the crystalline phase (with low resistance and high reflectance) and the amorphous phase (with high resistance and low reflectance) of PCM materials. The data writing and erasing operations are realized by the reversible phase transitions between the two phases. Since the operating stimuli also can be electronic or optical pulses, it can be used for both electronic and optical memory technologies.[25] PCM is firstly commercialized as optical memory including CD, DVD, and BD-RE.[6] These optical disks have been the main media for data storage for about 20 years. During the last two decades, the electronic devices, including computers, mobile phones, and so on, developed much faster than optical devices. Accordingly, the electronic PCM devices gradually became the main stream of PCM technology. Recently, however, the optically controlled phase-change memory restores the vitality due to the emergence of new applications, such as integrated all-optical non-volatile memory,[7] optically controlled in-memory computing,[8] and brain-inspired computing,[9,10] (see Fig. 1). The integration of optical memory on a chip [see Fig. 1(a)] demonstrates a practicable device architecture towards new types of memory or computing devices that beyond traditional optical disks. The in-memory computing [see Fig. 1(b)] means the recording and computing are realized in the same physical module, and thus can significantly enhance the computing efficiency and reduce the energy cost by avoiding the continual data transfer between the memory and the processor. In other words, it can break through the von Neumann bottleneck. The brain-inspired computing [see Fig. 1(c)] is mainly realized by fabricating the artificial neuros/synapses which are suitable to achieve the functionality of neuromorphic computing towards artificial intelligence and deep learning.

Fig. 1. The applications of phase-change memory in optical memory and in-memory/neuromorphic computing. (a) The integrated all-optical non-volatile memory based on PCM material Ge2Sb2Te5 (GST). Here, the PCM material is placed on the top of a waveguide. The data recording is realized by the optical near-field effects. The writing and erasing operations are achieved by the laser pulse passes through the waveguide while the read operation relies on the transmission signal modulated by the optical absorption of PCM materials. (b) The illustration of in-memory arithmetic computing by optically-controlled phase change. (c) The all-optical neuromorphic computing system based on PCM material GST. The input pulses can induce the phase change (recrystallization) of PCM material, which in turn changes the output signals. Both the in-memory computing and the neuromorphic computing utilize the accumulation effect of the crystallization. Panel (a) is reproduced from Ref. [7] with permission from the Nature Publishing Group. Panel (b) is reproduced from Ref. [8] with permission from the John Wiley and Sons. Panel (c) is reproduced from Ref. [9] with permission from the Nature Publishing Group.

Although the data processing at present is mainly dealt by electronic computers, the data transmission is dominated by optical communications in optical fibers. Therefore, repeated conversions between the electronic and the optical information are needed. If the data can be directly dealt by optical computing technology, the processing efficiency will be significantly enhanced. Moreover, the power consumption of optical computing also is believed to be lower than that of electronic computing.[11] As demonstrated by previous progresses,[79] optically controlled PCM should be a promising candidate for the advanced memory or computing technologies in the future. In addition, the laser-induced phase transition of PCM materials also shows potential in some other novel applications, such as the thermal emitter,[12] the thermal camouflage,[13] the plasmonic beam switching of metasurface and the bifocal lensing.[14] New opportunities are often companied by new challenges. Either in-memory computing or brain-inspired computing needs a fine control of the structures, properties and the phase-change dynamics of PCM materials. In this article, we review the understanding of laser-induced phase transitions in PCM materials which is the key for this technology. The dynamic process, the speed, the property change, the stability, the reversibility, and the cyclability of the transition closely influence the performances of PCM devices. Optical PCM involves both thermally and non-thermally driven phase transitions while electronic PCM mainly drives by Joule heating. Owing to the absence of violent thermal melting, the non-thermal phase transitions show some superior performances on the fine control of PCM materials by achieving unconventional phase transitions. However, the underlying mechanisms are still not fully understood due to the complicated light–matter interactions. Here, we mainly focus on the dynamic mechanism of laser-induced non-thermal phase transitions at atomic scale to provide fundamental supporting for optically controlled phase-change memory or computing.

2. Solid-to-solid amorphization induced by optical excitation

For optical PCM, the stimuli is laser pulse and the signal contrast (0/1) comes from the large difference in optical reflectance between the two phases (see Fig. 2).[15,16] Traditional understanding of the effect of laser pulse is heating.[5] During the amorphization process, a strong but short laser pulse heats the material over its melting point and thereby turns the material into its liquid phase. After the laser pulse irradiation, the liquid is rapidly cooled down to freeze the disordered atomic arrangement and thus transforms into its amorphous state. This process is called the melt-quenching process. For the recrystallization process, a weak but long laser pulse heats the materials over its glass-transition temperature and then triggers the spontaneous recrystallization process. The thermally driven phase transitions under laser irradiation explains the work principle of PCM quite well before.[2,5,1724] Thus, it has been believed to be the undoubted mechanism of PCM for several decades.

Fig. 2. The traditional mechanism of the melt-quenching and the recrystallization processes. The crystalline phase possesses a high reflectance while the amorphous phase possesses a low reflectance.

However, in 2011, by first-principles calculations, Li et al. proposed that a non-thermal effect of the optical excitation played important roles in the ultrafast laser-induced amorphization in PCM materials.[25,26] When electrons are excited from valence band to conduction band, the chemical bonds are weakened due to the electron transfer from bonding states to anti-bonding states. Such an instability mechanism also has been suggested in the femtosecond (fs) laser-induced disordering of Si and GaAs.[27,28] The results of the ab initio molecular dynamics (AIMD) simulation in Figs. 3(a)3(f) demonstrates that upon an excitation with an intensity of 9% valence electrons, the prototypical phase-change materials Ge2Sb2Te5 (GST) amorphizes within 24 picosecond (ps) at 700 K without thermal melting. In other words, the non-thermal effect of electronic excitation enables a solid-to-solid amorphization. This can also be confirmed by the mean square distance (MSD) during the AIMD simulations in Fig. 3(g), where the atomic diffusion of the solid-to-solid transition is negligible compared with that of the normal melting. Without violent thermal diffusions, the speed and the cyclability of the phase transition should be enhanced indicating better performances of PCM. AIMD simulations with different excitation intensities and temperatures also demonstrate that the amorphization temperature is in proportion to the excitation intensity. As displayed in Fig. 3(h), the higher the excitation intensity, the lower the amorphization temperature. At atomic scale, the excitation in phase-change materials shows an element-selective feature which makes Ge atom play a critical role in the phase changes. From the density of states (DOS) and the projected electron density in Figs. 3(i)3(k), it can be inferred that Ge atoms contribute more states than Sb atoms for the band-edge excitations. The excitation makes Ge atoms favor the sp3 hybridization instead of the p-bonding in crystalline phase.[26] Therefore, the excitation drives Ge atoms from the original octahedral sties to the new tetrahedral sites and thus triggers the amorphization in the way like “umbrella flip”.[18]

Fig. 3. Electronic-excitation-induced solid-to-solid amorphization in GST. (a)–(f) Structural snapshots during the AIMD simulation of amorphization under a 9% excitation. The crystalline structure turns to a disordered phase within 24 picoseconds (ps). (g) Mean square distances (MSD) during the AIMD simulations of the excitation-induced amorphization and the normal melting. The blue line indicates the exited state (Exc.) while the red line represents the melting state (Mel.). (h) The change of amorphization temperatures with different excitation intensities. The data are obtained from a series of AIMD simulations with different excitation intensities and temperatures. The red line is a quadratic fitting. (i) Partial density of states of GST. Upon a laser excitation, electrons are mainly excited from valence band maximum, i.e., the region labelled by “j”. Panels (j) and (k) show the projected electronic density distribution of the energy regions indicated by “j” and “k” in panel (i), respectively. Obviously, for the band-edge excitation, Ge atoms contribute more electronic states than Sb atoms. Panels (a)–(g) and (i)–(k) are reproduced from Ref. [26] with permission from the American Physical Society. Panel (h) is reproduced from Ref. [25] with permission from the John Wiley and Sons.

The non-thermal amorphization is also supported by some recent experiments. For example, Konishi et al. have reported that the fs laser irradiation in Ge10Sb2Te13 can result in an ultrafast (∼500 fs) decrease of the reflectivity [see Fig. 4(a)].[29] With relatively low laser fluences (12.8 mJ/cm2), the decreased reflectivity can recover gradually. When the fluence reaches 18.1 mJ/cm2, a long-lived decrease of the reflectivity can be observed indicating a permanent amorphization. Moreover, the scanning electron microscopy (SEM) measurement in Fig. 4(b) shows that the grain boundary contrast inherent to a polycrystalline structure remains unchanged. Therefore, they claimed that this is a non-thermal amorphization without melting. Siegel et al. also analyzed the time-resolved reflectivity of GST irradiated by fs laser [see Fig. 4(c)].[30] After irradiated by a relatively high laser fluence (60 mJ/cm2), the reflectivity rapidly decreases and then slightly recovers to the value of amorphous phase. Considering the fact that the reflectivity of the liquid-phase GST is lower than that of the amorphous GST,[31] the change of the reflectivity under the 60-mJ/cm2 irradiation is in consistent with the melt-quenching amorphization process. As such the reflectivity minimum should represent the value of liquid GST. In contrast, with a laser fluence of 36 mJ/cm2, the reflectivity directly turns to the amorphous value without reaching the liquid value, which implies a non-thermal amorphization without melting [see Fig. 4(c)]. In addition, compared with nanosecond (ns) laser, fs laser can induce the amorphization of GST with a faster speed and a lower fluence.[30] This phenomenon can also be explained by the non-thermal effect. For an ns laser, the pulse duration is much longer than the characteristic time of energy transfer from excited electrons to phonons. Therefore, thermal melting should be the dominant process. For an fs laser, the energy of excited electrons cannot be completely transferred to lattice within the pulse duration. Then, the above-mentioned non-thermal effect should play important roles in the amorphization of PCM materials with a temperature lower than the melting point.

Fig. 4. Femtosecond (fs) laser-induced amorphization of PCM materials. (a) Time-resolved reflectivity of Ge10Sb2Te13 after a 140-fs laser irradiation. The reflectivity intensity of PCM material decreases within 500 fs and then maintains for at least several tens of ps. (b) Scanning electron micrograph of a GST surface after irradiated by an 18.1-mJ/cm2 fs pulse. The dotted circle indicates the irradiated area where grain boundaries can be clearly seen. (c) Time-resolved reflectivity of the GST upon an fs laser irradiation. The data of different fluences are adopted from different areas of a light spot. The red dashed line indicates the reflectivity of liquid phase. (d) Time-resolved white line intensities of XANES of GST during the 600-ps laser-induced amorphization (red dots and line). Obviously, the intensities during amorphization are still far from the intensity of liquid state. Panels (a) and (b) are reprinted with permission from Ref. [29] ©The Optical Society. Panel (c) is reproduced from Ref. [30] with permission from the AIP publishing. Panel (d) is reproduced from Ref. [33] with permission from the American Physical Society.

Since optical reflectivity only indicates an electronic structure,[32] Fons et al. have used the white line intensity of x-ray absorption near edge structure (XANES) to detect the amorphization process via structural information.[33] Obviously, the results in Fig. 4(d) demonstrates that a sub-nanosecond laser makes GST change to its amorphous phase without turning into its liquid phase. The authors also attributed the phenomenon to a non-thermal mechanism. They suggested that the laser selectively excites the weaker bonds in GST which leads to the rupture of resonant bonding and thus results in the collapse of crystalline structure. Such a mechanism has also been proposed in GeTe by first-principles calculations, where a certain lattice distortion destroyed the resonant bonding and thereby spontaneously triggered a non-thermal amorphization.[34]

Recent years, time-resolved x-ray/electron diffractions have been utilized to directly catch the dynamic structural transitions. Waldecker et al. compared the time-resolved signals of dielectric function and x-ray diffraction (XRD) of GST under an fs laser irradiation [see Figs. 5(a) and 5(b)].[35] While the dielectric function decreases sharply after the pump laser pulse, the XRD intensity decreases gradually indicating a time-domain separation of electronic state from structural evolutions. The authors also estimated the lattice temperature via a Debye–Waller-type analysis based on an assumption that the phonons are close to a thermal distribution at all delay times. At the threshold laser fluence for amorphization, the estimated lattice temperature is close to the melting point of GST. According to these results, they suggested that the abrupt decrease of dielectric function comes from the ultrafast destruction of resonant bonding, while the structural amorphization is a result of thermal melting. However, this issue still needs further investigation, because the assumption of thermal phonon distribution and the above-mentioned non-thermal effect seem to be controversial. Although a certain degree of thermal effect is inevitable during a laser irradiation, whether the thermal melting can be completed within 5 ps is also questionable. Mitrofanov et al. also investigated the phase change dynamics of fs laser-irradiated GST by time-resolved XRD and x-ray absorption fine structure (XAFS).[36] A decrease of diffraction intensity without thermally-induced peak shift within 4 ps has been observed with indication of a non-thermal disorder effect. They proposed that the disruption of resonant bonding triggers local atomic distortion and thus leads to a long-lived intermediate state [see Fig. 5(c)]. The difference of XAFS signals between the intermediate state and liquid phase excludes the thermal melting process in this study. Hada et al. also studied the 400-nm fs laser-induced phase change dynamics in GST by time-resolved electron diffraction.[37] They suggested that the optical excitation-induces non-thermal amorphization of GST by triggering the motion of Ge atoms from octahedral to tetrahedral sites which is in consistent with the theoretical analysis.[26]

Fig. 5. Time-resolved (a) imaginary part of the dielectric function and (b) XRD intensity of GST after the fs laser irradiation. The intensity of the dielectric function dramatically decreases and then recovers to a certain degree which seems to be a behavior of electrons. In contrast, the intensity of XRD decreases slowly but remains stable which is regarded as the change of structures. (c) Schematic illustrations of the excitation-induced intermediate phase by structural rotation. Panels (a) and (b) are reproduced from Ref. [35] with permission from the Nature Publishing Group. Panel (c) is reproduced from Ref. [36] with permission from the Nature Publishing Group under a Creative Commons Attribution 4.0 International License (CC BY 4.0).
3. Mechanical effects of optical excitation

In addition to the local atomic motions, the optical excitation can also lead to global volume changes. For example, by time-resolved x-ray diffraction, Fons et al. detected a large volume expansion in the fs laser-irradiated Ge2Sb2Te5.[38] Their analyses suggested that the volume expansion cannot be explained solely by a thermal expansion. The first-principles calculations also proposed a giant expansion effect of electronic excitations. The excitation-induced expansion can be explained by the quantum electronic stress effect [see Figs. 6(a) and 6(b)].[39] By moving electrons from the valence band to the conduction band with atomic positions unchanged, the inner stress instantly appears due to the change of energy surfaces by the excitation. The global stress can be written as follows:[39]

Here, is the excited electron density, and is the effective deformation potential induced by an excitation. In GST, when the excitation intensity is larger than 6%, the effect of the stress is to drive the volume expansion [see Fig. 6(a)]. According to Eq. (1), the stress by excitation should be in proportion to the effective deformation potential and excitation intensity which has been confirmed by the calculations in a series of semiconductors [see Fig. 6(b)]. Since diamond and BN have very large deformation potentials, their excitation stresses are extremely large. Therefore, the excitation stress should be a general effect that is adapted for other materials.

Fig. 6. (a) The stresses and corresponding volume changes as a function of excitation intensity in rock-salt GST. Here, the negative stress leads to lattice expansions. (b) The stress as a function of the excitation intensity in a series of semiconductors. The inset shows the schematic illustration of inter-band excitation. CB and VB mean the conduction band and valence band, respectively. (c) The projected forces on individual atoms under a 15% excitation in rock-salt GST. The length of arrows indicates the magnitudes. (d) The magnitude of forces on Te atoms with different neighboring vacancies and excitation intensities in rock-salt GST. All panels are reproduced from Ref. [39] by permission of the PCCP Owner Societies.

To understand the rule of local atomic motions upon excitation, forces on individual atoms were studied by first-principles calculations in crystalline GST.[39] The local forces in rock-salt GST under the excitation display random directions [see Fig. 6(c)]. As such, fs-laser pulses usually make crystalline GST disordering. The underlying reason should be closely related to the unique structure for rock-salt GST, which have randomly distributed cations (Ge/Sb/vacancy) in the sub-lattices. For the magnitude of the atomic force on Te, it not only increases with the excitation intensity but also increases with the number of neighboring cation vacancies [see Fig. 6(d)]. In brief, both the direction and magnitude of the excitation force on an individual atom can be adjusted by the local symmetry of a material. If a particular structure can be designed, such as a superlattice GST, the excitation forces may be controlled with desired directions. These rules of atomic forces under optical excitation may offer some routes to control phase transitions to lower the power consumptions of PCM.[19,40]

4. Order-to-order phase transition induced by optical excitation

Although non-thermal phase transition has been demonstrated in some cases, the final result is usually an amorphization that is similar to the case by the thermal melt-quenching process. Surprisingly, in the most recent reports, Hu et al. and Matsubara et al. have observed an fs laser-induced rhombohedral-to-cubic phase change in GeTe by time-resolved electron diffraction and XRD, respectively [see Figs. 7(a)7(c)].[41,42] Such an order-to-order phase transition not only is helpful to understand the rule of initial atomic motions upon excitation but also shows the ability of laser controlling phase change in atomic scale that beyond the conventional amorphization. However, the mechanism of this phenomenon still has intensive debates. Hu et al. suggested that the phase transition is realized by direct motions of Te atoms along [001] axis [see Fig. 7(c)] and a subsequent shear deformation.[41] Matsubara et al. proposed that the cubic signal is an average effect of Ge-atoms rattling around the center positions under the excitation [see Fig. 7(d)].[42] Meanwhile, a previous theoretical calculation has suggested that the long and short bonds in GeTe will be randomly distributed by the excitations [see Fig. 7(e)].[43] This average effect of the random distribution could result in a cubic-like symmetry or a pseudo cubic phase. Therefore, the observed signal of rhombohedral-to-cubic transition has received quite different explanations.

Fig. 7. (a) Time-resolved XRD spectra of GeTe after fs-laser irradiations. The black lines are the XRD spectra before excitation for comparing. The merging of (003) and (101) peaks indicate the symmetry transition from rhombohedral to cubic. (b) The Debye–Waller factor deduced from the intensity of electron diffraction as a function of laser fluence. By comparing with theoretical simulations, the change of intensities is also demonstrated to be an indication of the symmetry transition from rhombohedral to cubic. (c) The schematic illustration of Te-atom displacement model where the cubic symmetry is a direct result of Te-atom movement. (d) The schematic illustration of Ge-atom rattling model where the Ge atoms dynamically move around the center of cubic sites. The average effect leads to a signal of cubic symmetry. (e) The schematic illustration of randomly distributed long and short bonds. The diffraction signal of such a structure also shows a cubic symmetry. Panels (a) and (d) are reproduced from Ref. [42] with permission from the American Physical Society. Panels (b) and (c) are reproduced from Ref. [41] with permission from the American Chemical Society. Panel (e) is reproduced from Ref. [43] with permission from the American Chemical Society.

To clarify the debates, Chen et al. performed first-principles calculations and real-time time-dependent density functional theory molecular dynamics (TDDFT-MD) simulations for GeTe by optical excitation.[44] TDDFT-MD is an advanced method to properly describe the excited-state behaviors because it contains the non-adiabatic electron–electron and electron–phonon couplings that beyond the traditional molecular dynamics simulations based on Born–Oppenheimer approximation.[45,46] The simulation with a 5% excitation intensity successfully reproduces the experimentally observed rhombohedral-to-cubic phase transition (Fig. 8). Due to the Peierls distortion, rhombohedral GeTe possesses a distortion along [001] axis that leads to three long bonds and three short bonds around a Ge or Te atom. When electrons are excited from valence band to conduction band, the electrons are mainly transferred from Te atom to Ge atoms [see Fig. 8(a)]. Then, the potential energy surface can be changed from a double-well shape to a single-well shape [see Fig. 8(b)]. As such, the increased gradient of potential energy instantly produces forces on Ge and Te atoms along the [001] axis or the z direction [see Fig. 8(d)]. Driven by the directional force, Ge and Te atoms coherently move to the opposite directions along [001] axis which activates the phonon modes. As a result, the long bonds and short bonds flip back and forth as shown in Fig. 8(e). Upon damping of the mode, the long and short bond lengths gradually transform to an equivalent value indicating a phase transition from rhombohedral to cubic phase. This is a truly rhombohedral-to-cubic transition rather than a pseudo one due to some sort of averaging effect. The temperature throughout the simulation is lower than the Curie temperature (700 K) of GeTe and thereby excludes the thermally induced phase transition [see Fig. 8(f)]. The atomic trajectory in Fig. 8(c) reveals that, the transition relies on collective motions without the traditional nucleation-and-growth process or atomic diffusion. Therefore, the phase-transition speed can be as fast as 1 ps.

Fig. 8. (a) The electron density distribution upon the 5% excitation in rhombohedral GeTe. The red arrows indicate the force acting on atoms. (b) The change of potential energy surface of GeTe due to electronic excitation. (c) The trajectories of atomic motions during the TDDFT-MD simulation of the phase transition in GeTe with a 5% excitation. Time evolution of (d) average force, (e) bond length, (f) temperature, (g) kinetic energy, and (h) momentum during the TDDFT-MD in GeTe under a 5% excitation. All panels are reproduced from Ref. [44] with permission from the American Physical Society.

Moreover, the kinetic energy and momentum of Ge and Te atoms that projected in three dimensions (x, y, z) show a strikingly higher component along [001] axis or z direction due to the excitation [see Figs. 8(g) and 8(h)]. This phenomenon clearly demonstrates the process is not the random thermal motion, which should have similar situations for the three different directions. In fact, the nominal temperature in Fig. 8(f) is a result of atomic motions including all the directions. Here, the authors recalculated the true thermal temperature just including atomic motions along x and y directions, which is not affected by the excitation [see the blue line in Fig. 8(f)]. The thermal temperature during the phase transition is clearly lower than 500 K. Therefore, the excitation-induced rhombohedral-to-cubic phase transition in GeTe is indeed a non-thermal process. Such a kind of order-to-order transition induced by fs laser not only represents a new type of phase transition but also supplies opportunities for exploring ultrafast phase change memory or some new optoelectronic applications.

5. Conclusion and perspectives

In this work, some recent progress on ultrafast laser-induced unconventional phase transitions in phase-change-memory materials are summarized to the community. Traditional understanding on the laser-induced amorphization is a thermally driven melt-quenching process. As the laser pulse duration becomes very short (e.g., femtosecond scale), the non-thermal effect of electronic excitation would play significant roles in the phase changes. First-principles calculations reveal the solid-to-solid amorphization in GST under optically electronic excitation. Analysis on electronic structures suggests that electronic excitation drives Ge atoms to move from the octahedral to the tetrahedral sites and thus triggers the amorphization very fast. Non-thermal amorphization of GST without melting is also supported by a series of experiments. Time-resolved reflectivity obtained by the pump-probe technology shows an amorphization process of GST without experiencing a liquid phase. Time-resolved XRD also shows locally structural changes of GST without melting, including local distortions and volume expansions. Furthermore, two recent experiments with time-resolved XRD and electron diffraction techniques both observed an order-to-order (rhombohedral-to-cubic) phase transition in GeTe by fs laser which is in contrast to the usual amorphization of PCM materials induced by laser. Due to the absence of atomic pictures in real time and real space, the underlying mechanism has been intensely debated. By TDDFT-MD, the atomic driving force and phase change dynamics are revealed. Owing to the change of potential energy surface upon electronic excitation, directional forces along [001] axis are instantly produced which directly drives the phase transition. The molecular dynamics simulations suggest that the order-to-order phase transition in GeTe can happen within 1 ps because it is completed by coherent atomic motions without the nucleation-and-growth process and the atomic diffusions. These unconventional phase transitions in PCM materials show the ability of ultrafast laser to manipulate structure beyond thermal melting. Compared with the conventional amorphization via a melt-quenching process, the solid-to-solid amorphization involves much less atomic diffusions. Moreover, the order-to-order phase transition not only possesses a speed with about 3 orders of magnitude faster than conventional phase-change random access memory, but also can control the property of phase-change materials with minimal atomic movements. Therefore, these unconventional phase transitions could be a promising route to promote optical storage and future optical computing with faster speed, lower power consumption, and better endurance.

Reference
[1] Ovshinsky S 1968 Phys. Rev. Lett. 21 1450 https://doi.org/10.1103/PhysRevLett.21.1450
[2] Breitwisch M J 2009 Phase Change Materials Raoux S Wuttig M Boston Springer 381 https://doi.org/10.1007/978-0-387-84874-7_17
[3] Sun H J Hou L S Wu Y Q Tang X D 2009 Chin. Phys. Lett. 26 024203 https://doi.org/10.1088/0256-307X/26/2/024203
[4] Sun H J Hou L S Wu Y Q Zhai F X 2008 Chin. Phys. Lett. 25 2915 https://doi.org/10.1088/0256-307X/25/8/050
[5] Wuttig M Yamada N 2007 Nat. Mater. 6 824 https://doi.org/10.1038/nmat2009
[6] Yamada N 2012 Phys. Status Solidi B 249 1837 https://doi.org/10.1002/pssb.201200618
[7] Ríos C Stegmaier M Hosseini P Wang D Scherer T Wright C D Bhaskaran H Pernice W H P 2015 Nat. Photon. 9 725 https://doi.org/10.1038/nphoton.2015.182
[8] Wright C D Liu Y Kohary K I Aziz M M Hicken R J 2011 Adv. Mater. 23 3408 https://doi.org/10.1002/adma.201101060
[9] Feldmann J Youngblood N Wright C D Bhaskaran H Pernice W H P 2019 Nature 569 208 https://doi.org/10.1038/s41586-019-1157-8
[10] Zhang Q Yu H Barbiero M Wang B Gu M 2019 Light: Sci. Appl. 8 42 https://doi.org/10.1038/s41377-019-0151-0
[11] Burr G W 2019 Nature 569 199 https://doi.org/10.1038/d41586-019-01406-0
[12] Du K K Li Q Lyu Y B Ding J C Lu Y Cheng Z Y Qiu M 2017 Light: Sci. Appl. 6 e16194 https://doi.org/10.1038/lsa.2016.194
[13] Qu Y Li Q Cai L Pan M Ghosh P Du K Qiu M 2018 Light: Sci. Appl. 7 26 https://doi.org/10.1038/s41377-018-0038-5
[14] Yin X Steinle T Huang L Taubner T Wuttig M Zentgraf T Giessen H 2017 Light: Sci. Appl. 6 e17016 https://doi.org/10.1038/lsa.2017.16
[15] Chen N K Li X B Wang X P Tian W Q Zhang S Sun H B 2018 Acta Mater. 143 102 https://doi.org/10.1016/j.actamat.2017.10.013
[16] Shportko K Kremers S Woda M Lencer D Robertson J Wuttig M 2008 Nat. Mater. 7 653 https://doi.org/10.1038/nmat2226
[17] Chen N K Li X B Wang X P Xia M J Xie S Y Wang H Y Song Z Zhang S Sun H B 2015 Acta Mater. 90 88 https://doi.org/10.1016/j.actamat.2015.02.015
[18] Kolobov A V Fons P Frenkel A I Ankudinov A L Tominaga J Uruga T 2004 Nat. Mater. 3 703 https://doi.org/10.1038/nmat1215
[19] Li X B Chen N K Wang X P Sun H B 2018 Adv. Funct. Mater. 28 1803380 https://doi.org/10.1002/adfm.201803380
[20] Akola J Jones R O 2007 Phys. Rev. B 76 235201 https://doi.org/10.1103/PhysRevB.76.235201
[21] Hegedus J Elliott S R 2008 Nat. Mater. 7 399 https://doi.org/10.1038/nmat2157
[22] Noé P Vallée C Hippert F Fillot F Raty J Y 2018 Semicond. Sci. Technol. 33 013002 https://doi.org/10.1088/1361-6641/aa7c25
[23] Zhu M Xia M Rao F Li X Wu L Ji X Lv S Song Z Feng S Sun H Zhang S 2014 Nat. Commun. 5 4086 https://doi.org/10.1038/ncomms5086
[24] Lin J C Long G Y Wang Y Wu Y Q 2013 Chin. Phys. B 22 124208 https://doi.org/10.1088/1674-1056/22/12/124208
[25] Li X B Liu X Q Han X D Zhang S B 2012 Phys. Status Solidi. B 249 1861 https://doi.org/10.1002/pssb.201200380
[26] Li X B Liu X Q Liu X Han D Zhang Z Han X D Sun H B Zhang S B 2011 Phys. Rev. Lett. 107 015501 https://doi.org/10.1103/PhysRevLett.107.015501
[27] Sokolowski-Tinten K Bialkowski J von der Linde D 1995 Phys. Rev. B 51 14186 https://doi.org/10.1103/PhysRevB.51.14186
[28] Stampfli P Bennemann K H 1990 Phys. Rev. B 42 7163 https://doi.org/10.1103/PhysRevB.42.7163
[29] Konishi M Santo H Hongo Y Tajima K Hosoi M Saiki T 2010 Appl. Opt. 49 3470 https://doi.org/10.1364/AO.49.003470
[30] Siegel J Gawelda W Puerto D Dorronsoro C Solis J Afonso C N de Sande J C G Bez R Pirovano A Wiemer C 2008 J. Appl. Phys. 103 023516 https://doi.org/10.1063/1.2836788
[31] Kieu K Narumi K Mansuripur M 2006 Appl. Opt. 45 7826 https://doi.org/10.1364/AO.45.007826
[32] Sun H B Wang L Wang H Y 2015 Physics 44 349 (in Chinese) https://doi.org/10.7693/wl20150601
[33] Fons P Osawa H Kolobov A V Fukaya T Suzuki M Uruga T Kawamura N Tanida H Tominaga J 2010 Phys. Rev. B 82 041203 https://doi.org/10.1103/PhysRevB.82.041203
[34] Kolobov A V Krbal M Fons P Tominaga J Uruga T 2011 Nat. Chem. 3 311 https://doi.org/10.1038/nchem.1007
[35] Waldecker L Miller T A Rude M Bertoni R Osmond J Pruneri V Simpson R E Ernstorfer R Wall S 2015 Nat. Mater. 14 991 https://doi.org/10.1038/nmat4359
[36] Mitrofanov K V Fons P Makino K Terashima R Shimada T Kolobov A V Tominaga J Bragaglia V Giussani A Calarco R Riechert H Sato T Katayama T Ogawa K Togashi T Yabashi M Wall S Brewe D Hase M 2016 Sci. Rep. 6 20633 https://doi.org/10.1038/srep20633
[37] Hada M Oba W Kuwahara M Katayama I Saiki T Takeda J Nakamura K G 2015 Sci. Rep. 5 13530 https://doi.org/10.1038/srep13530
[38] Fons P Rodenbach P Mitrofanov K V Kolobov A V Tominaga J Shayduk R Giussani A Calarco R Hanke M Riechert H Simpson R E Hase M 2014 Phys. Rev. B 90 094305 https://doi.org/10.1103/PhysRevB.90.094305
[39] Chen N K Han D Li X B Liu F Bang J Wang X P Chen Q D Wang H Y Zhang S Sun H B 2017 Phys. Chem. Chem. Phys. 19 24735 https://doi.org/10.1039/C7CP03103C
[40] Chen N K Li X B Wang X P Xie S Y Tian W Q Zhang S Sun H B 2018 IEEE T. Nanotechnol. 17 140 https://doi.org/10.1109/TNANO.2017.2779579
[41] Hu J Vanacore G M Yang Z Miao X Zewail A H 2015 ACS Nano 9 6728 https://doi.org/10.1021/acsnano.5b01965
[42] Matsubara E Okada S Ichitsubo T Kawaguchi T Hirata A Guan P F Tokuda K Tanimura K Matsunaga T Chen M W Yamada N 2016 Phys. Rev. Lett. 117 135501 https://doi.org/10.1103/PhysRevLett.117.135501
[43] Kolobov A V Fons P Tominaga J Hase M 2014 J. Phys. Chem. C 118 10248 https://doi.org/10.1021/jp412412j
[44] Chen N K Li X B Bang J Wang X P Han D West D Zhang S Sun H B 2018 Phys. Rev. Lett. 120 185701 https://doi.org/10.1103/PhysRevLett.120.185701
[45] Bang J Sun Y Y Liu X Q Gao F Zhang S B 2016 Phys. Rev. Lett. 117 126402 https://doi.org/10.1103/PhysRevLett.117.126402
[46] Runge E Gross E K U 1984 Phys. Rev. Lett. 52 997 https://doi.org/10.1103/PhysRevLett.52.997